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The Plasmodium of Physarum polycephalum is often used in the implementation of 
non-linear computation to solve optimization problems, and this organismal feature was 
not used in this analysis to compute perception and/or sensation in humans. In this paper, 
we focused on the Kanizsa illusion, which is a well-known visual illusion resulting from the 
differentiation-integration of the visual field, and compared the illusion with the adaptive 
network in the Plasmodium of P polycephalum. We demonstrated that the network pattern 
mimicking the Kanizsa illusion can be produced by an asynchronous automata-fashioned 
model of the foraging slime mold and by the real Plasmodia of P polycephalum. Because 
the protoplasm of the Plasmodium is transported depending on both local and global 
computation, it may contain differentiation-integration processes. In this sense, we can 
extend the idea of perception and computation. 
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1. INTRODUCTION 

A true slime mold, Physarum polycephalum (order Physarales, 
class Myxomycetes, subclass Myxogastromycetidae) has been used 
as a model organism to investigate biologically motivated uncon- 
ventional computing because the plasmodium of Physarum can 
solve various theoretical graph problems, including a Steiner tree. 
Spanning tree (Adamatzky, 2007), maze (Nakagaki et al, 2000; 
Nakagaki, 2001), and Voronoi diagram (Shirakawa et al., 2009; 
Shirakawa and Gunji, 2010). The plasmodium can also perform 
logical computing (Tsuda et al., 2004; Adamatzky, 2010) and cre- 
ate an optimal adaptive network (Tero et al, 2010; Gunji et al, 
2011). 

These computations, which produce particular patterns in 
relation to the input stimulus, were previously distinguished from 
psychological reactions such as sensation and emotion; how- 
ever, sensation and emotion are also considered higher- order 
computations in which pattern formation is applied to the devi- 
ation and/or representation of the input stimulus (Jones and 
Saeed, 2007; Adamatzky et al., 2013; Sakiyama and Gunji, 2014). 
In this sense, it is possible to say that even an amoeboid ani- 
mal can feel something (Humphrey, 2006). Whereas the retreat 
of the marginal front of the amoeba away from a repellant stim- 
ulus can be considered a lower-order computation, the reaction 
to a particular transportation of protoplasm from one margin to 
the other areas can be regarded as a sensation or a higher- order 
computation relative to the lower-order one. 

When the plasmodium of Physarum forages in an environ- 
ment, the extension or retreat of a local area can lead to the 
elongation and/or shortening of the tubular network and vice 
versa. This phenomenon reveals that a lower-order computa- 
tion can lead to a higher- order computation and vice versa. 
The synthesis of local and global information propagation can 



result in a reaction to an external stimulus coupled with a 
higher- order computation, which might be considered "percep- 
tion." This extended perception gives rise to the idea that the 
Plasmodium can emulate and/or compute the perception of a 
human. 

The visual perception of humans is not a simple local reac- 
tion to an optical stimulus, such as with a camera. Computing an 
optical stimulus is not only a local computation but also a global 
computation, and the local reaction to the optical stimulus can be 
modified depending on the global propagation of information. 
A visual illusion is one of the most intriguing examples of visual 
perception that results from the synthesis of local and global com- 
putation. In the Miiller-Lyer illusion, the apparent length of a 
local pattern can be shortened or extended depending on its sur- 
rounding pattern (Miiller-Lyer, 1896; Howe and Purves, 2005). In 
the aforementioned sense, biologically motivated computing can 
reveal visual illusion in terms of the analogy at the phenomeno- 
logical and mechanistic levels. We previously demonstrated that 
the distribution of foraging ants attracted to the honey dew along 
a Miiller-Lyer diagram can reveal a Miiller-Lyer illusion, wherein 
the ants in an attractive field can be compared to the neurons in 
a visual field (Sakiyama and Gunji, 2014). Not only ants but also 
the Plasmodium of Physarum can emulate various visual illusions 
because the plasmodium integrates lower- and higher- order com- 
putations. In this study, we focus on a particular visual illusion, 
Kanizsa illusory contours. 

Kanizsa illusory contours are subjective shapes induced from 
local pac-man elements (Kanizsa, 1979). According to Gestalt 
psychology, human vision is sensitive not to the input stimu- 
lus per se but to the relative values of input signals (Koffka, 
1935). This idea, when theoretically enforced (Ross and Pessoa, 
2000), explains that the Kanizsa illusion is triggered by the relative 
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contrast in the Kanizsa diagram, developed by contour comple- 
tion and finally completed by filling in (Pessoa et al, 1998). 

Although the process of contour completion was first 
explained using a simple collinearity-based model (Grossberg, 
1994; Grossberg and Pessoa, 1998), the disappearance of illusion 
cannot be explained by this model when the pac-man elements in 
the Kanizsa figure are replaced by crosses. The simple filling-in is 
unable to explain the perceived brightness and/or occlusion of the 
central area of the Kanizsa illusory contours, as noted by Coren 
(1972). To resolve these problems, a differentiation-integration 
approach is employed by some models (Gillam and Nakayama, 
2002). The edge of the bright area is detected by differentiation, 
and the contrast of the pac-man is also preserved in a memory. 
Through integration, this information creates a line connecting 
the pac-men, the brighter surface in the central area and the rel- 
ative depth featuring the occlusion cues (Kogo et al., 2010). By 
introducing a particular mechanism in which the local elonga- 
tions of contours interact and are canceled, the disappearance of 
the illusion for a particular variation of a Kanizsa figure can be 
explained. 

When the plasmodium of the Physarum is considered, the 
Kanizsa illusory contours might be emulated by a real forag- 
ing Plasmodium. We previously proposed an asynchronously 
updated automata- fashioned model for the plasmodium of the 
Physarum (Gunji et al, 2008b,a, 2011; Niizato et al, 2010). 
This model can mimic amoebic motions, adaptive networks and 
growth patterns in unstimulated environments of the plasmod- 
ium of Physarum. Because this model is based on the local 
detection of the environmental stimulus and the asynchronous 
transportation of protoplasm, the model can embed the synthesis 
of the differentiation-integration process. In this sense, the model 
and/or the real plasmodium of Physarum can connect separated 
pac-men by a tubular structure, which reveals an illusory contour. 
We ignore the brighter surface and occlusion cues of the central 
illusory area in the Kanizsa illusion and focus only on the con- 
tour creation. Because the filling-in and the relative depth can 
be simulated if our Physarum model is coupled with a particular 
mechanism, the significance of our model and the real Physarum 
should not be lost. 

In this paper, we demonstrate that our model for foraging 
Physarum can reveal Kanizsa illusory contours if the pac-man 
faces the center, whereas the illusory contours disappear if the 
pac-man faces outward. In this result regarding stimulation, 
the pac-man is implemented as an attractant in the automata- 
fashioned model. We also demonstrate that the same patterns 
are obtained in the experiment involving the real plasmodium of 
the Physarum. In this case, the pac-man is implemented using a 
nutrient-rich dry agar sheet. These observations suggest that the 
transfer of biological information resulting from the synthesis of 
local and global interactions can create what is called extended 
visual perception, thus allowing visual illusion, i.e., biological 
and/or unconventional computing can emulate the perception of 
humans. 

2. MATERIALS AND METHODS 
2.1. DEFINITION OF THE MODEL 

Models of Physarum plasmodium have been suggested. For 
example, weighted graph model (Tero et al., 2006, 2007), particles 



model (Jones, 2010), and network growth model (Takamatsu 
et al, 2009), and so on. We use the Vacant- Particle (VP) model 
proposed by Gunji et al. (2008a,b) to explain the adaptive net- 
work of the Plasmodium of Physarum. This is an asynchronous 
cellular automaton defined on a discrete lattice plane in which 
the asynchronous application of transition rules provides amoe- 
boid motility. This model has the advantage of not only sim- 
ulating the amoeboid cellular motion of Physarum but also 
forming adaptive networks that depend on the surrounding envi- 
ronment (Niizato et al., 2010; Gunji et al., 2011; Tani et al., 
2011). 

First, we explain the basic algorithm of the model. The addi- 
tional rules coupled with the model are defined later. The coordi- 
nates of a site on the discrete lattice plane are described by a pair 
of integers (i, j). For all sites (i, j), a binary value c(z, j) is assigned. 
c(i, j) = 1 indicates that the site (i, j) is filled with the plasmod- 
ium and implies that the site is inside of the plasmodium. These 
are called internal sites. Similarly, c(i, j) = 0 indicates that the site 
(i, j) is outside the plasmodium. These are called external sites. In 
particular, an internal site that confronts at least one external site 
in a von Neumann neighborhood is defined as a membrane. In 
the left column of Figure 1 the normal internal sites are yellow, 
and the membrane sites are orange, whereas the external sites are 
not colored. 

We use the symbol N to indicate an aggregation of internal 
sites. 

N={{i,j) \c{i,j) = l} (1) 

Similarly, the symbol M indicates an aggregation of membrane 
sites. 

M = {{ij) eN\c{i + lj) = Ov c{i - 1,;) = Ovc{iJ + l) 
= Ovc{iJ - l)= 0} (2) 

Next, we explain the processes that provide amoeboid cell motility 
for N. Note that the word "neighborhood" means "von Neumann 
neighborhood." 

The temporal transition of a plasmodium from time t to time 
t + 1 is illustrated in the left column of Figure 1 . There are four 
child processes in the background of this change, and these are 
depicted in the right part of Figure l(i-iv). 

(i) We randomly choose one site from M. This selected site is 
called the stimulus point (SP). The SP is colored blue in 
Figure l(i). 

(ii) Because the SP is a membrane cell, there is at least one 
external site in the neighborhood of the SP. We also ran- 
domly select an external site from among these neigh- 
bors. This selected external site is called the growth point 
(GP). The protoplasm at SP then moves to GP. This move- 
ment means that the value c at the GP changes to 1. 
Instead of this movement, an external void inserts into 
the SP. The void that has entered the plasmodium is 
called the VP 

(iii) The VP migrates between the internal sites N by exchang- 
ing places with neighboring internal sites. This sequential 
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FIGURE 1 I Basic algorithm of the model. The left colunnn indicates the tinne evolution of the model fronn tinne f to f + 1 . The right-hand portion indicates the four 
child processes: (i) the choice of stinnulus point; (ii) the growth to the outside and entry into the VP; (iii) the nnigration of VP with a gel transition; and (iv) the escape of VP 



movement of the VP mimics a protoplasmic flow. The 
most important condition is that the protoplasm experi- 
encing an exchange itself produces a gel transition. Because 
of this, these cells are never exchanged with the VP again. 
Therefore, the migration pathway of the VP never intersects 
with itself The red-colored sites in Figure l(iii) indicate the 
gelled protoplasm. The VP continues searching exchange- 
able cells in the neighborhood and exchanging places if 
possible. 

(iv) In the case where the VP has no exchangeable neighborhood 
sites or the VP has exchanged a particular given num- 
ber of times, the VP escapes from the plasmodium. This 
site is called the escape point (EP). Then, the value c at 
the EP becomes 0. Eventually, we observe that the proto- 
plasm that was at the EP has moved to the GP. Therefore, 
these processes conserve the number of internal sites, i.e., 
the covered area by plasmodium does not change tempo- 
rally. Finally, the entire gelled protoplasm undergoes a sol 
transition. 



random selection of the SPs and GPs prohibits unidirec- 
tional growth. Hence, to acquire real Physarum plasmodium- 
like behavior, some taxis must be defined regarding this 
model. In fact, the introduction of the localization of prob- 
abilities to the random choice process yields taxis for the 
model. 

2.2. ADDITIONAL RULES 

It is well known that real Physarum plasmodium has 
a negative taxis to white light. Thus, we propose an 
Avoid Bright region VP (ABVP) model as defined 
below. 

For all sites (i, j) on the lattice plane, a new binary value b(i, j) 
is defined. This value h{i,]) corresponds to the brightness of the 
site. h(i,]) = 1 indicates that the site is bright, and b(i,j) = 0 
indicates that the site is dark. 

In the basic model, there are three types of random choice 
processes, which are detailed below. 



This is the basic algorithm of the model. These pro- 
cesses are called the "one time step," i.e., they occur all at 
once. The characteristic of this model is a gel transition. 
This condition maintains the integrity of the plasmodium 
and enables amoeboid motion. However, this basic model 
Physarum cannot move over a long distance because the 



1 . The choice of the SP. 

2. The choice of the GP. 

3. The choice of the direction of movement of the VP. 

For each of these processes, we introduce variations of 
probabilities. 
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2.2. 1. The choice of the SP 

First, for all sites (i, j), we define the integer value s(i,j) depend- 
ing on c(i,j) and b(i,j). The value s(i,j) is determined as foUows 
[Figure 2(i)]: 



0 (c(/,i) = 1) 

1 {c{i,j)=OAb{i,j) = l) 

2 {c{i,j)=OAb{i,j)=0). 



Next, the new value w(i, j) is defined for all of the membrane sites 
depending on s(i, j), as defined in the equation below. 

w{i,j) = s{i + 1, j) + s{i - 1,;) + s{ij + l) + s{i,j - l) . (4) 

Here, w(i,j) reflects the local environment of the site This 
w(i, j) value determines the probability p5(i, j) that the membrane 
site (i, j) will be selected as the SP. 



Ps{ij) 



w(m, n) 

(m,n)eM 



(5) 



Therefore, the membrane sites that have many external and dark 
neighborhoods are selected with a high probability. 

2.2.2. The choice of the GP 

Let a membrane site (i, ;) be selected as the SP. The GP is then 
determined using the probabilities described below. 

For a neighbor site of the SP (i, j), which for convenience we 
write as the probability that the site will be 

selected as the GP is determined by the value s that was previously 
defined in [Figure 2(ii)]. 



s{i + l,i)+s{i - l,j)+s{i,i + l) + s{i,j - l)' 

(6) 

Note that if a neighboring site is internal, then the probability 
choice becomes 0 because of the definition of s. 

2.2.3. The choice of the direction of movement of the VP 

The above two conditions promote biased growth toward the dark 
sites. Next, we introduce another rule to implement repellency 



from the bright region. To decrease the plasmodium under the 
bright area, the VP is expected to move toward the brighter area. 
Let the VP be at the site (i, j); then for all neighbor sites of 
the VP, a value r(i^ /) is defined as follows [Figure 2(iii)]: 



0 (^0"'/) = 0V (/",/) is gelled) 

(3) r(/",/) = ] 1 {c{i\f) = 1 A is not gelled A b{i\f) = O) . 

2 {c{i\f) = 1 A {i\f) is not gelled A b{i\f) = l) . 

(7) 

The probability pni(i\ f) that the neighbor site (i^ /) will be cho- 
sen as the next destination of the VP is then determined as in the 
following equation: 



r{i + 1, j) + r{i - + r{i,j + l) + r(i, j - l) " 

(8) 

Note that if the value r(^^/) = 0 for all neighborhood 
values, i.e., the probabilities pm(i\f) = 0 for the entire neigh- 
borhood, then the escape process occurs according to the basic 
algorithm. 

However, the effects of this condition are very restrictive in the 
sense that only local information of the VP is used, and this cor- 
rection is valid only when the VP is at the edge of a bright/ dark 
region. 

2.2.4. Linear structure correction 

An internal site (i, j) is called a "linear structure" if both sides of 
(i, j), (i + 1, j), (i — 1, j) or + 1), 1) are external sites. 

During the process of randomly choosing a stimulus point, if 
a candidate site (i, j) is a linear structure, then a positive constant 
L is added to w(i,j). From previous work (Gunji et al, 2011), it 
is known that a larger value of L produces a straighter linear net- 
work. However, a value for L that is too large will invalidate the 
negative taxis for the bright regions. Therefore, in this paper, the 
value of L was set to 10. 

2.2.5. Thicl(ness rules 

In the basic algorithm, the number of internal sites does not 
change throughout the development. This implies that the area 
covered by the model plasmodium is unchanged. A real Physarum 
Plasmodium, however, changes its cell thickness spatiotemporally 
to spread widely or concentrate in a narrow area. In particular. 




Frontiers in Cellular and Infection Microbiology 



www.frontiersin.org 



February 2014 | Volume 4 | Article 10 | 4 



Tani et al. 



Kanizsa illusion in slinne nnold 



under the condition that we consider here, a real plasmodium 
would likely converge upon dark regions, and its cell thickness 
would increase. 

To implement this behavior, we add the next two rules to our 
model. Let (ig,jg) be the GP and (ie,je) be the EP at the time 
step t. 

After the escape process, 

1. ifb(ig,jg) = 0 and b(ie,je) = 1? then c(ig,jg) = 0, 

2. if h(ig,jg) = 1 and h(ie,je) = 0? then c(ie,je) = 1- 

What do these rules mean? (1) If the protoplasmic flow moves 
from a bright region to a dark region, then protoplasm is used to 
increase the thickness at the SP. Accordingly, the growth process 
appears to be halted [Figure 3(i)]. (2) In contrast, if the proto- 
plasm moves from a dark region to a bright region, then the cell 
thickness at the dark region decreases instead of protoplasmic 
spreading at the bright region. This rule invalidates the escape 
process [Figure 3 (ii)]. 

The former rule decreases the area covered by the Plasmod- 
ium, and the latter rule increases this area. These rules are applied 
deterministically depending on only the brightness at the GP and 
the ER 

2.2.6. Fluctuation 

The final additional rule is related to fluctuations. At each time 
step, the definition of the value s, which was given in section 2.2. 1, 
is changed with the probability ]y such that 



0 (c(/,i) = l) 

1 (c(/,i) = 0) 



(withpy) . 



(9) 



The linear structure correction is not applied here. Therefore, the 
stimulus point is chosen depending on only the local information 




regarding the state of the sites. Because the SP is chosen uniformly, 
this rule accordingly promotes growth in the bright region, which 
results in the reconstruction of the networks. As a result, shorter 
networks are given. In this paper, we set pf = 0.05. 

2.3. EXPERIMENTAL CONFIGURATION 
2.3.1. The model 

A detailed description of our Physarum model is provided herein. 
Pac-men in the Kanizsa diagram are defined as dark regions, and 
other areas are defined as bright regions (Figure 4A). For the con- 
trol experiment, the pac-men in the Kanizsa diagram are turned 
180° (Figure 4B). Clearly, in the case of the control, the illusory 
contour does not appear in the perception of humans. The resolu- 
tion of both of these images is 120 x 120 pixels. The lattice plane 
on which the models are defined has the same size. 

We conducted experiments 100 times for each image and 
examined the network topology obtained after five million time 
steps for each experiment. The maximum limit for the number of 
VP movements for each time step was set at 300. The initial state 
of the Plasmodium was determined using the following method: 

1. For all sites on the lattice plane, we defined c = I; 

2. A site that has a value c = I was randomly chosen as the SP; 

3. The VP enters the plasmodium at SP, although growth does 
not occur; 

4. The VP migrates between the internal sites conforming to the 
rules of the VP model and eventually stops; 

5. The site at which the VP stops becomes c = 0 (Because the 
growth process was halted, the area covered by the plasmod- 
ium decreases by one site); 

6. Repeat 2-5 until the entire covered area is reduced to <45% of 
the area of the entire lattice plane. 

This method provides cancellate networks that cover the entire 
plane with many medium-size holes. These networks imitate the 
real Physarum tubular networks that are formed on an unstimu- 
lated agar plate. 

In addition to simulation studies, we conducted a real plas- 
modium experiment in which pac-men in the Kanizsa diagram 
are implemented with nutrient rich agar (attractant). In the next 



FIGURE 3 I Thickness rules, (i) If GP is dark and EP is bright, then 
c(GP) = 0, and the covered area decreases. Conversely, (ii) if GP is bright 
and EP is dark, then c(EP) = 1 , and the covered area increases. 
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Kanizsa diagram 




Control 


FIGURE 4 1 Stimulus figures are given as geometrical patterns of 
brightness (black or white). (A) Kanizsa diagrann and (B) control figure. 
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section, we describe the experimental configuration for a real 
Plasmodia experiment. 

2.3.2. Real Plasmodia 

We cultivated Physarum plasmodia using the method of Camp 
(1936). Briefly, the plasmodia were cultured on wet paper towels 
laid out on glass Petri dishes in a plastic box. The space below the 
towels was filled with tap water to maintain humidity. Oatmeal 
flakes were provided daily as food. We maintained the culture 
environment in the dark at 25°C. 

In this study, we introduced a sheeted food source made from 
nutrient-rich agar to make attractant figures. The sheet was easily 
shaped into various figures and sizes. The nutrient-rich sheet was 
produced as follows: 400 ml of a solution containing 1.5% agar 
and 20 mg/ml oatmeal powder was poured into a 295 x 220 mm 
plastic tub; the resulting sheet was dried thoroughly in a natural 
convection oven at 50° C. 

In the experiment, we made three pac-men from the above 
nutrient agar sheet. The radius of a pac-man was 5 mm, and 
the distance between the centers of two pac-men was 15 mm 
(Figure 5). We used a 1.5% unamended agar plate in a 9 cm diam- 
eter Petri dish as an experimental field. We put the agar sheet 
pac-men onto the plain agar plate in a fixed arrangement and 
kept it at rest until the sheets were rehydrated (approximately 
1 h). Next, we inoculated 5 mg of plasmodial fragments onto the 
center of each of pac-man and photographed the plasmodia at 
intervals of 30 min. We captured images using a single-lens dig- 
ital camera (EOS Kiss X6, Canon, Japan) from directly above, 
and the experimental system was kept at 25° C. We repeated the 
experiments 24 times and performed the same number of control 
experiments. 

3. RESULTS 

3.1. MODEL RESULTS 

We classified the resulting networks based on their topology. The 
networks were divided into four types (Shirakawa and Gunji, 
2007) (Figure 6) as follows: 

• Type I is a triangle that connects three pac-men; 

• Type II is a V- shape obtained from a triangle lacking of side; 

• Type III has another node that links with all of the pac-men 
and one side of a triangle; 

• Type IV is a simple star structure. 



A B 




FIGURE 5 I Experiment using real Physarum (scale bar: 10 mm). The 

radius of a pac-man was 5 mm, and the distance between the centers of 
two pac-men was 15mm. (A) The captured image. (B) Obtained by 
subtracting the background of (A). (C) Adjusted the brightness. (D) 



However, for the Kanizsa diagram, there were a few 
undifferentiated networks that could not be classified any 
of the above types. Therefore, all of networks these were cat- 
egorized together as "others" and excluded from the following 
statistical test. Obviously, type I networks complement the 
illusory contour of a Kanizsa diagram. 

The frequencies of each type are given in Table 1. Clearly, in 
the case of the Kanizsa diagram, type I networks were obtained 
more frequently than the control case. Type IV networks were 
not obtained at all, while this type comprised over 40% of the 
resulting networks in the control case. A x -square test indicated 
that there was a significant difference between these two data 
[X-square(3, N = 191) = 56.7, p < 0.001]. Therefore, this dif- 
ference in the frequencies between the Kanizsa triangle and the 
control figure reflects the difference in the configuration of the 
pac-men. 

What was the cause of this difference? To clarify this, we illus- 
trate an example of the time development of the network for each 
stimulus image (Figures 7A,B). In the case of the Kanizsa dia- 
gram, the illusory contour is very stable. It formed completely 
before t = 2, 000, 000 and was maintained until t = 5, 000, 000. 
In contrast, in the control, the network changed over time and was 
not stable. This is a general disposition of the types in Table 1. 
According to this perspective, the explanation of the difference 
is a difference in the stability of the Physarum network for each 
stimulus image. This suggests the presence of a unique potential 
structure for each image. The lack of a mouth for the pac-man 
may work as some sort of "energy barrier" for the plasmodial 
networks. In the Kanizsa diagram, these barriers produced the 
triangular networks because all of the pac-men were arranged 
face-to-face. However, the configuration of the control figure did 
not work similarly. In this sense, the resulting networks reflected 
the global configuration of the pac-men. 

3.2. RESULTS USING REAL PLASMODIA 

In these experiments, we compared real Physarum networks with 
our model. Because the real plasmodia networks usually have 
more complex structures than our model (Figure 8), the classi- 
fication of network types is not obvious, unlike the model cases. 
Therefore, we adopted another analytical method to estimate the 
spatial distribution of the plasmodia. 

The captured images of the real plasmodia were processed 
using Image J image processing software (Rasband and Image J, 



C D 




Binarization. All processes were executed by using the ImageJ software. 
In the binary image, the plasmodia are indicated by black pixels. The 
threshold of binarization is determined not to delete delicate structure of 
Plasmodium. 
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FIGURE 6 I Classification of the resulting networks according to their topology. Type I is a triangle; type II is a V-shape; type III is a slingshot-like shape; 
and type IV is a star or Y-shape. A type IV network was not obtained for the Kanizsa diagrann. 



Table 1 | The frequency of network types for two stimulus patterns. 



Stimulus image 


Type 1 


Type II 


Type III 


Type IV 


Others 


Total 


Kanizsa diagrann 


79 


9 


3 


0 


9 


100 


control 


46 


3 


8 


43 


0 


100 



1997-2012). Figures shows each step of the binarization pro- 
cess for the image. Figure 5A is the original photograph, and 

(B) is obtained using "Subtract Background", a macro function 
of Image J, on (A). The "Adjust/Brightness and Contrast" com- 
mand to automatically normalize the brightness produces image 

(C) . The threshold of binarization is determined according to the 
brightness value of the site on which the plasmodia are thinly 
spread to maintain the delicate structures. Obvious noise sur- 
rounding the Plasmodium, such as dust or bubbles in the agar 
plate are removed manually. Finally, the binarized image (D) is 
obtained. The time series of plasmodial networks are shown in 
Figures 7C,D, and some examples of the binarized images of the 
networks are shown in Figure 8. According to these images, the 
Plasmodium rarely entered the central area of the illusory triangle 
in the case of the Kanizsa diagram. In contrast, the corresponding 
areas of the control figures are occupied by the plasmodium. 

To clarify this, we produce Figure 9 by stacking binarized 
images for all of the trials. Because the initial configuration has 
rotational symmetry, the images that are 120° and 240° rotations 
of the original are also used. In Figure 9, the brightness of each site 
indicates the intensity of lapping, which is the probability of the 
Plasmodium being at the site. The sites that are occupied more 
frequently have higher brightness levels. Therefore, the whole 
image represents the distribution of the probabilities. We have 
used the "Lookup table/Yellow Hot" color chart in ImageJ. To 
improve the visibility of the figures, a filter was applied to cut 
off the weak signals. This threshold is determined as 10% of the 



theoretical maximum of the signal, and this value is common to 
both conditions. To remove the effects of Kanizsa illusory contour, 
the pac-men are covered with black in the lower row. Obviously, 
the type I -like structure is apparent in the Kanizsa diagram case 
(Figure 9A), and the structure that is the sum of type I and type 
IV networks appears in the control (Figure 9B). This observation 
does not conflict with our model results because type I networks 
appeared dominantly in the Kanizsa diagram case, but in the 
control, type I and type IV networks were obtained in almost the 
same ratio. Next, we emphasize this difference numerically. 

We defined the target region of evaluation as depicted in 
Figure 10. The yellow- colored marginal area of the fan that cen- 
ters a pac-man was divided into six regions. All of these fans had 
the same central angle of 15° and were labeled A to F accord- 
ing to the direction from the center. The same labeled areas 
were identified. For each region, the ratio of the area occupied 
by Plasmodium was measured. We precisely counted the num- 
ber of black pixels, which indicated that the site was occupied 
by Plasmodia in the binary converted images (see, for example. 
Figure 5D). To normalize the values, the number of black pixels 
was divided by the total number of pixels for each region. The 
relationships between the mean values of the occupation ratios 
and directions are plotted in Figure 11. 

In the Kanizsa diagram cases, the plasmodia converged to areas 
A and F. In contrast, in the control cases, the Physarum spread 
uniformly, and there was no clear peak. In the model case, there 
were significant differences in the occupation ratios between the 
two figures for all areas at a 0.1% level according to the Wilcoxon 
rank sum test (A: W = 74, 996, p < 0.001; B: W = 34, 774, p < 
0.001; C:W = 26, 542,p < 0.001; D: W = 26, 952,p < 0.001; E: 
W = 37, 559, p < 0.001; F: W = 74, 480, p < 0.001). Similarly, 
in the case of actual plasmodia, we could see that there were 
significant differences at the same level for only areas C and D, 
which occupied the center of fan (A: W = 2795.5, p = 0.416; 
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FIGURE 7 I Typical evolution of the networks of model (A,B) and real 
Plasmodium (C,D) (scale bar: 10 mm). (A) For the Kanizsa diagrann, the 
network structures are highly stable. (B) In the control diagrann, the network 
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changes gradually and beconnes shorter. (C) The plasnnodiunn does not enter 
the central area of the illusory triangle. (D) The plasnnodia occupy the nniddle 
area of the pac-nnen. 



B: W = 2402, p = 0.448; C: W = 937, p < 0.001; D: W = 1085, 
p < 0.001; E:W = 2335, p = 0.304; V: W = 2753, p = 0.52). In 
this regard, the results of the model and the real plasmodia are 
highly consistent. Namely, the illusory contour of the Kanizsa dia- 
gram is complemented by the plasmodia in both cases. These 
results are derived from the directed cell motility depending on 
the configuration of the attractant. 

4. DISCUSSION 

These results demonstrate that our visual perception of the illu- 
sory contours of the Kanizsa diagram can be emulated by both the 
asynchronous automata- fashioned model and by real plasmodia 
as an analogy at the phenomenological level. This finding does 
not mean that the plasmodia have human-like intelligence or per- 
ceptive faculties. It is clear that Physarum plasmodia cannot be a 
direct model of our brain because plasmodia do not have even 



primitive nervous systems. Nevertheless, we can say that there is 
abstract homogeneity between plasmodial behavior and nervous 
activity in our brain. This analogy has meaning at a mechanistic 
level. 

It is known that our consciousness contains perception and 
sensitivity arising from neural information that is processed in 
our brains. The causality between particular activities of neu- 
rons and our consciousness has not been sufficiently investigated. 
Nevertheless, we can regard this nervous activity as a type of 
"computation" that is calculated by complex neuronal firing pat- 
terns. In other words, we can say that our consciousness is a 
product of a program consisting of multiple neural networks. 
Obviously, we cannot consider that the program is the same as an 
ordinary computer program at the syntactical level. However, the 
semantics of the program can be discussed. Taking this perspec- 
tive, it is possible that the perception arising from brain activity is 
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FIGURE 8 1 The resulting networks from the real Physarum 
experiment (scale bar: 10 mm) corresponding to the Kanizsa 
diagram (A) and the control (B). The tinne stannp of each 



networks are shown in the innage. Note that these networks are 
the results of individual trials and these photos are not in a 
tinne series. 






FIGURE 10 I Target region of analysis. The marginal areas of the fans that 
are colored yellow are divided into six regions and nanned A-F according to 
the direction. For each divided region, we counted the number of black 
pixels and estimated the area covered by the Plasmodium. An occupation 
ratio of a region is given as the ratio of covered area to the total area of the 
region. 



Kanizsa diagram 



control 



FIGURE 9 I The upper diagrams are produced by stacking 
binaralized images of the plasmodial networks for all trials. (A) 

Kanizsa diagram and (B) control. The brightness of each site 
indicates the intensity of lapping; that is the existence probability of 
the Plasmodium at the site. A filter that cuts-off the weak signals 
was applied to improve visibility. This threshold is common to both 
figures. In the lower diagrams, to wipe out the effect of Kanizsa 
illusory contour, the pac-men are covered by black (scale bar: 10 mm). 



analogous to plasmodial behavior, which is of a different material 
and complexity altogether. 

To perceive the illusory contours of the Kanizsa diagram, 
both the differentiation and integration processes are necessary. 
Depending on the local gradient of the stimuli, detection of the 
edge of the pac-man corresponds to a differential process, and the 
combination of the local data based on their global configuration 
corresponds to an integral one. Therefore, the former is a lower- 
order process that finds the difference in stimuli in a relatively 
narrow area on the retina, and the latter is a higher- order 
process that integrates this information and constructs a visual 
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FIGURE 11 I The relationship between the mean value of the occupation 
ratio and direction. For the Kanizsa diagram, there are two peaks at region A 
and F in both experinnents. In the nnodel case, type I and type III networks 



field. This means that the illusory contours appear with a sort of 
information contraction from lower to higher orders. 

What plasmodial behavior corresponds to the each of differ- 
ential and integral processes? Plasmodia search their surrounding 
environment and react to the attractant and the repellant. 
This is differential. For integration, plasmodia propagate local 
information using protoplasmic flow to decide their overall 
behavior. In our model, the imbalance of the choice of SP 
depending on the local state and the movement of VP based on 
asynchronous application of the transition rules corresponds to 
the differential and integral processes, respectively. Therefore, the 
model and real plasmodia have analogous functionality to the 
processes of visual perception of the illusory contours. In this 
way, plasmodia can combine lower and higher- order informa- 
tion smoothly without a specialized pivotal organ, such as a nerve 
center, that oversees the entire system. 

It is important to note that the resulting networks' comple- 
mentary illusory contour cannot be obtained as an optimized 
network at the global level. Actually, we can link all of the pac- 
men with more shorter paths. For example, the V-shaped network 
(type II) and the Steiner tree (star shape, type IV) have more total 
shorter paths than the triangular network (type I). However, these 
networks were not generally obtained in either experiment. If the 
plasmodia oversee the whole and choose optimized networks, tri- 
angular networks would not appear. This means that the illusory 
contours seem to result from the redundancy of the network, 
which is considered to be an error or mistake in the Physarum 
computer studies that consider some optimization problems. 

In the Kanizsa diagram, the mouth of a pac-man works as an 
obstacle to change shorter paths and provides metastable net- 
works such as in type I because all pac-men are face-to-face. 
However, in the control, all pac-men are back-to-back, and the 
mouths of the pac-men do not behave as barriers. Then, the 
transition of the network structure to a more stable and shorter 



canceled each other out. Therefore, a clear peak does not appear. There is 
also no peak in the real plasmodiunn case because of the isotropic growth of 
the plasnnodiunn. 



network occurs quickly. This indicates that the resulting networks 
reflect the potential structure of each spatial configuration of 
pac-men. 

The use of an analogy to recognize that plasmodial behav- 
ior and visual perception have the same functionality provides 
us with one perspective about the origin of illusory contours. 
When we pick out information from our visual field, it is impor- 
tant to determine the borders of the objects using gradients 
of stimuli such as brightness or color. However, it is ineffi- 
cient to process all local information at the same level because 
the volume of information is too huge. Therefore, to efficiently 
find the border, we must compress the local information to 
some extent. In this way, we can say that illusory contours are 
a by-product of the information contraction corresponding to 
the moderate optimization of the plasmodial networks in our 
analogy. 

We have equated the plasmodial behavior and our perception 
at the mechanistic level and discussed the analogous homo- 
geneity between them. Finally, we refer to a point of difference 
between our study and previous studies concerning Physarum 
computation. Whereas we used the true slime mold Plasmod- 
ium as a computational device for a visual perception model, 
in conventional works, plasmodia have been used for optimiza- 
tion problems. We can say that the latter viewpoint considers 
Physarum computation to be an effective bio-inspired algorithm, 
similar to a genetic algorithm or a neural network. 

Our perspective is completely different. We are attempting to 
understand the computational processes of organisms through a 
model based on bio-computation. A bio-computer is the most 
suitable device for modeling the computational processes of an 
organism. To understand biological computational processes, it is 
necessary and important to construct models of these processes 
using a bio -computer and to discover the abstract connections 
between the model and the phenomena. Our experiments are a 
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Starting point for Physarum computational studies that suggest 
the computabihty of our "perception" and "sensitivity." 
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